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CN ■ Abstract. A chiral quark model is described which is regularized in terms of Lorentz 

invariant non-local interactions. The model is regularized to all loop orders and it 
ensures the proper quantization of the baryon number. It sustains bound hedgehog 

Qv ' solitons which, after suitable centre of mass corrections, can adequately describe the 

f~>». . nucleon. 

o 

5^ : I SOME SPECIFICITIES OF CHIRAL QUARK 

5 : MODELS 

^, 

'p'' This work was done in collaboration with Wojciech Broniowski from Krakow. 

Oh. We consider chiral quark models which encompass three sectors. The vacuum and 

r-| ! soliton sectors, which are treated in the mean-field (leading order in A^^) approx- 

imation, and the meson sector, which describes the (next to leading order in N^ 
vibrations of the vacuum sector. Not all models are applicable to the three sectors. 
For example, constituent quark models, in which quarks interact with confining 

5^ \ forces, cannot describe the vacuum sector, that is, the Dirac sea. However, they 

can and do describe the excited states of baryons, a thing which the chiral quark 
models cannot do (except, possibly, the A) for lack of confinement. 

Chiral quark models (nor any of the other low energy quark models) have not 
been derived from QCD. The only serious attempt to derive them from QCD is the 
instanton gas model [1,2]. In this approach, the chiral quark model is derived by 
calculating the propagation of quarks in a gas of instantons. A regularized effective 
theory results, as it should. It predicts both the value of the cut-off and the form 
of the regulator. The non-local regularization discussed here has the same form as 
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the one derived from the instanton gas model. Unfortunately, the quark models 
derived from the instanton structure of the vacuum do not lead to quark or color 
confinement. This serious limitation serves as a reminder that we have not really 
succeeded in deriving low energy effective theories from QCD. 

Other so-called "derivations" of quark models from QCD involve more guesswork 
than derivation. Most telling is their inability to derive a regularized model. If 
infinities appear in an effective theory, one should seek the physical processes which 
prevent the infinities from occurring. Invoking the roughly 200 MeV QCD cut-off 
is not a serious argument. Nor does QCD imply in any sense that the quark-quark 
interaction at low energy should be a one-gluon exchange with a modified gluon 
propagator. The regularizations used so far in the Nambu Jona-Lasinio type models 
for example (proper-time regularization being the most commonly used one so far), 
are nothing but renormalization techniques in which a finite cut-off is maintained. 
Not only is this arbitrary but such regularizations are flawed with problems. 

One might argue that the value of the cut-off should not matter. Indeed it would 
not if the effective theory consisted, for example, in eliminating some high energy 
degrees of freedom and using the remaining degrees of freedom to work out the 
dynamics of low energy phenomena. In such a case, one might expect the cut-off to 
be much larger than the inverse size of the composite particles and the results not 
to be sensitive to the cut-off. In chiral quark models, however, this is not the case. 
The cut-offs required to fit /,r are about 700 MeV, hardly larger than the p or the 
nucleon mass. This is a fact of life, whether we like it or not. One can of course 
simply discard such models, but better models do not seem to be forthcoming. 

II THE SOLITON IN THE NON-LOCAL CHIRAL 

QUARK MODEL 

The non-local chiral quark model is defined by the euclidean action: 
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In this expression, F^ = (l,i'~^^Ta), q{x) = {x\q) is the quark field, and r is a 
regulator. The regulator is assumed to be diagonal in momentum space and it 
has a range which defines an effective euclidean cut-off A. For example, we could 

take (A;|r|A;') = Sk^k'f^ {k"^) with r [k"^) = e~2A^, where /c is a euclidean 4-vector 
kfj, = (uj,k] with k'^ = co'^ + 1? . The interaction term of the action (1) can be 
viewed as a contact 4-fermion interaction involving the delocalized quark fields: 

i) {x) = {x \r\ q) = d^y (x \r\ y) q (y) . (2) 

An action of the form (1) is derived from the instanton gas model of the QCD 
vacuum [1,2], which predicts a cut-off function of the form: 
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where p is the instanton size. The the cut-off is determined by the inverse instanton 



size p. The form (3) has r {z = 0) = 1 and r 



However, at large 



euchdean momenta k, the form (3) is no longer valid and the cut-off function is 
dominated by one gluon exchange. It decreases as p- (with possible logarithmic 
corrections) and not as p-. We find that the fall-off of the regulator at large 
euclidean k'^ does not affect the soliton properties very much. For this reason, we 
have felt free to use various simple forms of cut-off functions, such as a gaussian, 
which have an additional advantage in that they can be analytically (although 
arbitrarily) continued to negative values of k'^. We shall see below that the analytic 
continuation is required to include the valence orbit. Similar regularization has 
been used by the Manchester group [3] in the meson and vacuum sectors. Various 
regularization schemes are reviewed in chapter 6 of Ref. [4]. 

The euclidean action allows us to calculate the partition function Z = 
J D (a) D (a'^j e~ v"'" '' and the ground state energy E = —jglnZ. The partition 

function cannot be written in the form Z = Tr e~^^ because the regulator in the 
action (1) prevents us from defining a hamiltonian H. We are also unable to quan- 
tize the quark fields but we shall see that the baryon number is nonetheless properly 
quantized. 

We work with the equivalent bosonized form of the action: 
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in which case the partition function is given by the path integral Z = J D (cp) e ^^'^\ 
We refer to ipoTa = S + i'j^TaPa, as the "chiral field" and we say that the chiral 
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field is "on the chiral circle" if, for all x, we have S"^ 
is an x-independent constant mass. 

We have calculated a localized and time independent stationary point of the 
action (4), consisting of a chiral field with a hedgehog shape S (r) + i'-y^XafaP (r) 
[5]. The shape of the fields and the soliton energy can be calculated in terms of 
the energies ex {oo) of the quark orbit. The "Dirac hamiltonian" is diagonal in the 
energy representation, although it remains energy dependent. The quark orbits 
\uj, \uj) satisfy the equations : 
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The energy of the soliton is: 



Esol = NcC^ai + TT ^^^ Yl '■ — : — ^^1~T + KFPJ / ^3X >^l (x) - vac. (6) 

ZTT J-oo Y^ too + ex [uj) 2G^ J 

where —vac. means that we subtract the vacuum energy. In the vacuum, P = 0, 
5* = Mq and there is no valence orbit contribution e^ai- The latter is discussed in 
the next section. 



Ill THE QUANTIZATION OF THE BARYON 
NUMBER AND THE VALENCE ORBIT 

We calculate the baryon number from the Noether current associated to the 
gauge transformation q (x) -^ e~*"*^^^g (x). It turns out to be: 

B = -TT d^y. '". ^ • (7) 

2mNc J-oo j;^ tuj + ex (uj) 

The extra term '^^' in the numerator arises from the fact that the regulator r 
does not commute with a {x) . Its effect is to make the residues of all the poles of 
the quark propagator ^^ , ^ (^^ equal to unity. This effectively quantizes the baryon 
number in a manner which does not seem to be related to the topology of the 
hedgehog field. ^ This is most fortunate because, a priori, there is no reason to 
expect a theory, in which we cannot quantize the quark field, to yield a properly 
quantized baryon number. 

The expression (7) suggests a way to include the valence orbit so as to ensure 
that the baryon number of the soliton, relative to the vacuum, is equal to unity. 
We calculate "on-shell" pole of the quark propagator in the hedgehog background 
field by searching for a solution of the equation iu + ex (a;)|^^^g = 0. Because of 
the regulator, the solutions are scattered all over the complex u plane. However, it 
is well known that, in the local theory, where we set r = 1, and for a hedgehog field 
with winding number unity, a well separated bound orbit with grand spin and parity 
O"*" occurs with energy e^ai close to zero [6]. In the non-local theory, we find that a 
solution of the equation u = ie^ai iy^) can always be found on the imaginary oj axis, 
close to the origin u; = 0, and that no other pole occurs in the vicinity. We therefore 
ensure that the soliton has a baryon number i? = 1 by deforming the integration 
path over uj in such a way as to include the contribution of this pole. This requires 
an analytic continuation of the regulator. Such a continuation is arbitrary but the 
analytic continuation does not extend as far from the origin as e^ai- Indeed, since the 
soliton size is small, /c^ > is large and this, on the average, makes k'^ = — e^^^ + k"^ 
less negative. Unfortunately however, the form (3) of the regulator, predicted in 
the instanton model, does not allow any analytic continuation whatsoever, thereby, 
strictly, prohibiting its use in the soliton calculation. 



^) Nor is the soliton stabilized by the topology of the chiral field. 



IV RESULTS OF SELF-CONSISTENT SOLITON 

CALCULATIONS 

The model parameters are the couphng constant G appearing in the lagrangian, 
the cut-off A appearing in the regulator and the current quarlc mass m. The 
values of the three parameters are constrained by fitting the pion decay constant 
/^ = 93 MeV and the pion mass to m.^ = 139 MeV. The expression used to calculate 
the pion decay constant f-,, is: 



f^ = 2NfN,M, 
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valid in the chiral limit m — *> and it is not identical to the Pagels-Stokar formula 
[7]. This leaves one undetermined parameter which we choose to be the constituent 
quark mass Mq at zero 4-momentum. The pion decay constant /tt sets the scale. 
Grossly, soliton energies increase and soliton radii diminish as /tt increases (see 
table 2). 

Figure 1 shows the soliton energy Esoi as a function of the free parameter Mq. 
A soliton is a bound state of A'^c = 3 quarks which polarize the Dirac sea. With 
a gaussian regulator, it is formed if Mq ^ 276 MeV, that is, for a sufficiently 
strong coupling constant G ^ 4.7 x 10"'^ MeV~^. The bound state occurs when the 
energy of the system is lower than the energy N^Mq of N^. free constituent quarks 
in the vacuum: Egoi < NcMq. The mass Mq is the on-shell constituent quark mass, 
obtained by searching for the pole of the quark propagator in the vacuum. It is the 
solution of the equation k"^ + (r|Mo + m) =0, which requires an analytic 

continuation of the regulator to negative values of /c^. Figure 1 also shows N^Mg. 
At the critical value Mq ^ 276 MeV, the two curves merge. The contribution N^e^ai 
of the valence orbit is also shown. At the critical value of Mq, the energy e^ai of the 
valence orbit, which is the on-shell mass of a quark propagating in the hedgehog 
field, becomes a well distinguished bound orbit. 

At Mq ^ 309 MeV, the curve displaying NcMq on figure 1 abruptly stops. Indeed, 
for larger values of Mq, the poles of the quark propagator no longer occur for real 
values of k"^. This means that quarks can no longer materialize on-shell in the 
vacuum. This feature is discussed in chapter 6 of Ref. [4] and it has been considered 
by several authors as a sign of quark confinement [8-10]. In fact, when a pole of 
the quark propagator disappears from the real k"^ axis, it simply moves into the 
complex plane. Such poles indicate instability of the assumed vacuum state against 
the addition of a single quark. 

However, our calculation shows that, in the background soliton field, the on-shell 
valence orbit continues to exist and so does the soliton. Unfortunately, the regulator 
also introduces extra unwanted poles in the propagators of colorless mesons, so 
that the model does not express color confinement. Similar unwanted poles occur 
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FIGURE 1. The energy of the sohton [in MeV] (bold soUd hne), Nc times the free-space quark 
mass (sohd hne) and the valence contribution to the soliton energy (dashed line) plotted as 
functions of the parameter Mq [in MeV]. A Gaussian regulator is used; A (dots) is fitted to 
U = 93 MeV. 



in proper-time regularization [11]. Our ignorance as how to continue propagators 
in the complex k'^ plane reflects our ignorance of the conflning mechanism [12]. 

Apart from the solitons consisting of three valence quarks we flnd stable solitons 
consisting of a single valence quark in the background soliton fleld (see flgure 2) 
as well as of two valence quarks. Similar solutions have been found in the linear 
sigma model with valence quarks [13]. 

Figure 3 shows the scalar and pseudoscalar flelds S (x) /Mq and P (x) /Mq of 
the soliton obtained with several values of Mq, together with the soliton quark 
density p (x). Note that, within the soliton, the flelds do not lie on the chiral circle 



and S^ (x) + P2 



X 



<Ml 



Indeed, the pion component P (x) never reaches the 



values — Mq. This is a new dynamical result. This is the only calculation, as far 
we know, in which one can check dynamically whether the chiral fleld remains or 
not on the chiral circle. It could not be checked in the renormalized linear sigma 
model, because close lying Landau poles occur which make the soliton unstable 
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FIGURE 2. The energy per quark [in MeV] for the sohton with three valence quarks (bold 
line), the soliton with one valence quark (dashed line) and the free-space quark mass Mq plotted 
as functions of the parameter Mq [in MeV] . 



against high gradients in the fields [14,15]. It could also not be checked in local 
theories which use proper-time regularization because, in such theories, the soliton 
is unstable unless the fields are constrained to remain on the chiral circle [16,17]. 
No such instability occurs with the non-local regularization. 

The soliton we obtain with non-local regularization has a structure which lies 
midway between a Friedberg-Lee soliton [18,19] (in which the pion field has a 
vanishing classical value), and a Skyrmion [20,21] (in which the chiral field is con- 
strained to remain on the chiral circle). This raises the problem of the collective 
rotational motion of the soliton. If the deformation in spin and isospin space is 
stable enough to sustain a rotation without significant distortion, then the A can 
be described as a rotation of the soliton and the A^ — A mass splitting can be 
estimated by cranking. If, however, the deformation is small, the A may be better 
described as a bound state of quarks with aligned spins and isospins. We have not 
tackled this problem yet. 

Table 1 shows some properties of calculated solitons for various values of the 
mass parameter Mq. Rather good values of qa are obtained. The soliton mass and 
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TABLE 1. Properties of self-consistent soliton solutions obtained with a gaussian 
regulator. 

energies need to be corrected for spurious centre of mass motion (see table 2). 

The fields which describe the soliton break translational symmetry. The center 
of mass of the system is not at rest and it makes a spurious contribution both to 
the energy and to the mean square radius (more generally, to the form factor). 
This spurious contribution is not measured and it should be subtracted from the 
calculated values. The subtraction occurs at the next to leading order (in A*";,) 
approximation. A rough estimate can be obtained from an oscillator model. If N^ 
particles of mass m move in a Is state of a harmonic oscillator of frequency hu, 
the centre of mass of the system is also in a Is state and it contributes ^huj = 
(P^) /2Ncm to the energy. We have therefore corrected the soliton energies by 
subtracting (P^) /2Esoi from the calculated energy. Furthermore, in the oscillator 
model, the center of mass contributes a fraction -^ of the mean square radius, so 
that we have corrected the mean square radius by multiplying the calculated value 
by a factor equal to f 1 — -^ j . 

Table 2 shows the result. The soliton energies and radii are then considerably 
closer to the experimental values observed in the nucleon. 
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TABLE 2. Elimination of spurious cm. motion. 
Gaussian regulator, A fitted to /^ — 93 MeV; * A 
fitted to ,U = 1.25 X 93 MeV. 



V CONCLUSION: WHY TAKE THE TROUBLE? 



The non-local regularization effectively cuts out of the quark propagators the 4- 
momenta which are larger than the cut-off. The non-local regularization makes the 
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FIGURE 3. Self consistently determined fields and baryon densities {Aur^ p) for various values 
of Mq; a gaussian regulator is used. 

theory finite at all loop orders. The simpler proper-time and Pauli-Villars regular- 
ization schemes regularize the quark loop only and they require extra independent 
cut-offs when next to leading order meson loops are included. Both the real and 
the imaginary parts of the action are regularized, while the anomalous properties 
remain independent of the cut-off [22-24], and the baryon number remains prop- 
erly quantized. In proper time and Pauli-Villars regularization schemes only the 
real part of the action is regularized and the imaginary part is left unregularized 
in order to enforce correct anomalous processes. Why not limit the 3-momenta of 
the quarks, thereby avoiding unwanted extra poles in the propagators? Breaking 
Lorentz covariance in the meson sector is annoying in that it requires to boost 
composite particles calculated in their rest frame. 
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